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The practice of garden bird feeding is a global phenomenon, involving millions of 
people and vast quantities of food annually. Many people engage in the practice of 
feeding assuming that birds gain some benefit from the food they provide, yet recent 
studies have revealed the potential for detrimental impacts as well. However, there is 
still a paucity of information on the impacts of feeding, including the ubiquity of these 
impacts among and within feeder-visiting species. Consistency in feeder use among 
birds is likely an important determinant of this. Individual birds and species that make 
frequent use of feeders are more likely to experience both the benefits and detrimental 
impacts of supplementary food. We investigated patterns of feeder use by garden birds 
visiting experimental feeding stations in Auckland, New Zealand, with the specific aim 
of determining whether use of supplementary food was consistent or variable among 
individuals and species. We used camera traps as well as Radio Frequency Identification 
(RFID) technology to examine intra- and interspecific feeder visitation patterns and to 
discern species associations. Eleven bird species were detected using feeding stations, 
however, two introduced species (house sparrow Passer domesticus and spotted dove 
Streptopelia chinensis) dominated visitation events. These species were present at 
feeders most frequently, with the largest conspecific group sizes. Significant associations 
were detected among a number of species, suggesting interspecific interactions are 
important in determining feeder use. We also found within-species differences in 
feeder use for all focal species, with individual variation greatest in house sparrows. 
Furthermore, season had an important influence on most visitation parameters. The 
observed individual and species-specific differences in supplementary food resource use 
imply that the impacts of garden bird feeding are not universal. Crucially, particularly given 
the avifaunal context in New Zealand, resource dominance by introduced species could 
have potential negative outcomes for native species conservation in cities. 


Keywords: competition, individual variation, interspecific interactions, resource use, supplementary feeding, 
urban wildlife 


INTRODUCTION 


Garden bird feeding is a phenomenally popular activity in many parts of the world, including 
in New Zealand, Australia, the UK, Europe, and the USA (Jones, 2017). Participation rates for 
engagement in bird feeding are consistently estimated at between one- to two-thirds of households 
(Jones and Reynolds, 2008; Galbraith et al., 2014; Orros and Fellowes, 2015b). In recent decades bird 
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feeding has shifted from a predominantly winter-only activity 
to a pastime commonly practiced year-round (Jones and 
Reynolds, 2008; Horn and Johansen, 2013; Galbraith et al., 2014). 
Effectively a massive ecosystem-scale intervention, bird feeding 
has numerous potential implications for the biology and ecology 
of feeder-visiting birds (Jones, 2011), as well as the wider faunal 
community (e.g., Bonnington et al., 2014; Orros et al., 2015). 
Although studies of bird-feeding impacts in urban habitats are 
rare, there is mounting evidence to confirm that garden bird 
feeding can be profoundly influential for urban-dwelling bird 
communities (Amrhein, 2014). For example, feeding can alter 
body condition, reproductive outputs, adult survival, disease 
dynamics, community assemblages, and migration (Robb et al., 
2008; Jokimaki and Kaisanlahti-Jokimaki, 2012; Galbraith et al., 
2015, 2017; Orros and Fellowes, 2015a; Plummer et al., 2015; 
Wilcoxen et al., 2015). The vast body of scientific literature on 
the influence of supplementary feeding on a wide range of species 
and non-urban habitats corroborates many of these findings (e.g., 
Boutin, 1990; Clout et al., 2002; Ilarri et al., 2008; Schoech, 2009; 
Ruffino et al., 2014). 

Whether the impacts of supplementary feeding are universal 
among and within feeder-visiting species remains largely 
unstudied. Certainly, food availability is an important factor 
acting to limit bird populations (Newton, 1980), affecting 
reproductive success and survival of many bird species 
in different systems (Martin, 1987). Despite the additional 
deliberate (i.e., bird feeding) and unintentional (e.g., refuse) food 
resources available to birds in urban systems, demand for food 
can be high due to high bird densities reducing the per capita 
amount of food available (Seress and Liker, 2015). This demand is 
illustrated by an experimental study in Arizona, USA, that found 
supplementary food in urban areas was depleted much faster than 
equivalent amounts provided in natural habitats (Shochat et al., 
2004). Thus, competition for food resources in urban systems, 
including supplementary food, may be high. The competitive 
ability of animals, in discovering and dominating a resource, 
plays an important role in the structuring of numerous faunal 
assemblages (e.g., in ants; Parr and Gibb, 2012; Bertelsmeier et al., 
2015). Intraspecific and interspecific asymmetries in competitive 
abilities commonly give rise to dominance hierarchies (Holway, 
1999), dictating resource access and consequently determining 
which individuals or species gain benefits from the resource 
(French and Smith, 2005). Dominant competitors may displace 
others via interference (physical exclusion from a resource via 
aggression) or exploitative (rapid discovery and removal of 
a resource) competition (Bertelsmeier et al., 2015). Agonistic 
interactions are certainly common at bird feeders, with body 
size found to be a critical factor in determining outcomes 
(Tamm, 1985; Wojczulanis-Jakubas et al., 2015). Speed of novel 
food discovery (Tryjanowski et al., 2015) and bird densities at 
feeding locales (Galbraith et al., 2015) can also vary by species. 
Furthermore, in natural habitats, different foraging strategies 
are frequently used by different individuals within a species to 
optimize foraging efficiency (Gustafsson, 1988). Consequently, 
it is unlikely that the ability to exploit supplementary food 
is consistent among all feeder-visiting birds, both within and 
between species. 


Intraspecific and interspecific differences in use of 
supplementary food could have important implications for 
population- and community-level impacts of urban bird 
feeding, potentially determining the mechanisms by which 
changes to survival, reproduction, migratory patterns, and 
community organization occur (Newton, 1980; Robb et al., 
2008). Few studies, though, have looked at individual variation 
in supplementary food use or how species associations at feeders 
affect resource access in an urban context specifically (but see 
Cowie and Hinsley, 1988; Crates et al., 2016; Jack, 2016). Bird 
populations in urban habitats are subject to different pressures 
than those in more natural environments, frequently resulting 
in differences in ecology, behavior, and life history (Chace and 
Walsh, 2006; Chamberlain et al., 2009; Seress and Liker, 2015; 
Garcia et al., 2017; Lepczyk et al., 2017). Furthermore, urban 
areas, particularly in New Zealand, are hotspots for introduced 
(i.e. nonnative/exotic/alien) bird species (Day, 1995; Duncan 
et al., 2003; van Heezik et al., 2008; Spurr, 2012; Davis et al., 
2014), with bird feeding implicated in the success of some of 
these species (Strubbe and Matthysen, 2007; Peck et al., 2014; 
Orros and Fellowes, 2015b). Thus, it is immensely important 
to study supplementary food use in situ, in urban areas where 
most bird feeding occurs, to gain a realistic understanding of 
the demand for these resources, the urban-specific competitive 
interactions which may be occurring, and the wider implications 
for urban birds. 

Here we explore avian visitation patterns at experimental 
bird-feeding stations established in the gardens of volunteer 
households in Auckland, New Zealand. Via this experiment we 
have also investigated the impacts of bird feeding on avian disease 
dynamics (Galbraith et al., 2017) and avian community structure 
(Galbraith et al., 2015). We know from the latter that our feeding 
regime significantly altered bird communities at feeding locales, 
prompting a shift toward communities heavily dominated by 
introduced birds, primarily house sparrows (Passer domesticus) 
and spotted doves (Streptopelia chinensis). This study investigates 
whether there were intraspecific and interspecific asymmetries in 
feeder use that could be indicative of resource dominance, giving 
insight into the community-level observations of Galbraith et al. 
(2015). We also consider whether supplementary food use is 
modified by interspecific interactions or seasonality. The demand 
for supplementary food resources, and the associated competitive 
interactions at feeders, are likely to vary across seasons (e.g., 
Ottoni et al., 2009; Cox et al., 2016) due to fluctuations in natural 
food availability and physiological (thermoregulatory as well as 
reproductive) demands on birds. We used camera traps at feeding 
stations to identify species-level patterns of feeder use and 
examine species associations, and Radio Frequency Identification 
(hereafter “RFID”) technology to explore the feeder-visitation 
patterns of individuals. Specifically, the objectives of our study 
were to: (1) examine whether feeder use varies among and 
within species to determine whether birds exploit supplementary 
food equally; (2) determine whether individuals and species are 
consistent in their use of supplementary food over time or if 
use varies seasonally; and (3) explore the associations between 
species at feeders that may modify access to supplementary 
food. 
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METHODS 


Experimental Feeding Stations 

Experimental feeding stations were established at 11 urban 
residential properties in northern Auckland, New Zealand, as 
part of a wider study of typical bird feeding practices (Galbraith 
et al., 2015, 2017). The study area is largely suburban residential, 
with a population density of 1,600/km” in 2006 (New Zealand 
Census data, www.stats.govt.nz). Properties representative of the 
study area were selected from a pool of 42 volunteered properties. 
Feeding stations (a low feeding table (40 x 80 cm, 17-cm high), 
a seed feeder, and a mesh bread-tube; Supplementary Figure 
1A) were active for 18 months from March 2012 to September 
2013, and householders were responsible for provisioning them 
on a daily basis for the duration of the study. The feeding 
regime consisted of 4-5 slices of bread and 1 metric cup 
of birdseed (white millet, Hungarian millet, hulled oats, and 
canary seed blend) per day; householders were asked to put 
the food out between 0700 and 0800 h NZST. This regime 
and the design of the feeding stations reflected typical feeding 
practices of the New Zealand public (Galbraith et al., 2014). The 
majority (75%) of feeding participants in New Zealand throw 
food directly onto the ground, rather than using structures or 
containers (Galbraith et al., 2014); however, for the present study 
it was necessary to use a fixed structure and food containers 
to standardise experimental feeding among properties and to 
enable data collection. Properties were a minimum of 900 m apart 
(Supplementary Figure 2), with no detections of banded birds at 
properties other than their place of capture. For full details of the 
experimental setup see Galbraith et al. (2015). 


Species Visitation and Association Data 

The visitation patterns of feeder-visiting species and species 
associations were examined using camera traps. Camera 
trap data were collected at all feeding stations (n = 11) 
over four sampling periods: Austral winter 2012 (June-July; 
non-breeding season), spring 2012 (October-November; early 
breeding season), summer 2013 (February; late breeding season), 
and winter 2013 (June-July). In each sampling period, three 
cameras (ScoutGuard SG570V, HCO Outdoors, USA; 1.5 m from 
the feeding station, at a height of 0.5 m) were rotated around 
the feeding stations over a 4-week timeframe, and operated for 
4-6 nights at each feeding station depending on weather. Three 
full days of recordings were scored for analysis. The cameras 
were programmed to record upon motion-activation, with a 
10-min delay period following a recording event to increase 
independence of observations. Three consecutive photographs, 
stamped with the date and time, were taken at 2-s intervals 
upon motion activation. These photo sets were considered as one 
observation event. 

We used co-occurrence at feeding stations to identify 
associations among species. For scoring, we counted the 
maximum number of individuals of a species observed 
simultaneously in an event. For analysis we included only those 
recordings with at least one “feeder visitor,’ henceforth referred 
to as a “visitation event.” Birds on the feeding station structure 
itself (i.e. on the feeding table or feeders) were considered 
feeder visitors; birds on the ground or elsewhere nearby were 


disregarded. We also noted the level of food remaining at the time 
of the visitation event: bread, seed, or both food types remaining 
(visible in the containers) or none (no food visible, negligible 
quantities available). 


Individual Visitation Data 

To investigate individual patterns of feeder use we used 
RFID and Passive Integrated Transponder (hereafter “PIT”) 
tag technology. At three study properties birds were captured 
via mistnetting as part of our wider study of feeding impacts 
(for full details see Galbraith et al., 2017). We PIT-tagged 
a subset of captured individuals, and fitted RFID antennae 
to feeding stations (Supplementary Figures 1A, 2). The RFID 
reader setup (Microchips Australia Pty Ltd) consisted of a 
30 x 25 cm (inner dimensions) coil antenna attached via an 
RS232 serial cable to a LID-650N decoder (Trovan Ltd, UK) 
mounted on nearby posts that also supported a 20 W solar 
panel (Supplementary Figure 1B). RFID readers were powered 
with a 26A h, 12 V sealed lead acid battery (HAZE Solar 
Gel, USA) in weatherproof housing. Acrylic walls supported 
antennae and encouraged birds to pass through the antenna to 
reach the food containers (Supplementary Figure 1A). Antennae 
were intentionally designed to allow multiple birds access to the 
feeding stations simultaneously, reflecting a typical bird-feeding 
situation. Readers recorded the individual code (ID) of all PIT- 
tagged feeder-visitors and the time and date of the visit, with a 1-s 
read-delay between consecutive reads of the same ID, and were 
active between 0600 and 2100 h NZST from 22 September 2012 
until the end of the study (13 September 2013). Readers did not 
distinguish between arrival/departure movements. 

We PIT-tagged five species (Supplementary Table 1) under the 
conditions of our animal ethics permit (University of Auckland 
Animal Ethics Committee Permit R921). Birds were tagged over 
three capture rounds conducted at 6-monthly intervals, with the 
monitoring duration (number of days from PIT-tag implantation 
to the end of the experiment) among returning individuals 
varying between 134 and 377 days (mean = 291.5 d + 10.7 SE). 
Sample sizes depended largely on capture rates of those species. 
PIT-tagging took place on only those days where two experienced 
researchers were present in the field. Only adult birds, without 
apparent injuries or clinical signs of illness, were PIT-tagged. 
We used Trovan Unique ID100 implantable PIT tags (2.12 x 
11.5 mm; 0.1 g; Trovan Ltd, UK) with a unique 5-byte code. 
These were injected subcutaneously in the back of birds above 
the scapula following the methods of Nicolaus et al. (2008). The 
tags are pre-sterilized by the manufacturer, and come ready to 
use with a disposable needle. During each procedure, one person 
cleaned the insertion site and injected the PIT tag, while a second 
person held the bird and gently pulled up the skin at the top of 
the back to facilitate injection. The perforation of the skin at the 
insertion site was closed by applying a small quantity of surgical 
adhesive (Vetbond, 3M, St. Paul, MN, USA). This method has 
been used on small passerines without obvious negative effects on 
the birds (Nicolaus et al., 2008; Nomano et al., 2014). Five PIT- 
tagged individuals were recaptured during subsequent days or 
capture rounds and we found all tags to be in place with no visible 
problems for these individuals (Supplementary Figure 1C). 
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For analysis of PIT-tag data we omitted readings from the 
same individual that occurred within 2s of each other, to obtain 
a more conservative estimate of feeder visitation rates. Two key 
parameters were calculated for analyses: the total number of 
reads per day for each individual, and the presence/absence of 
each individual at the feeder for each day that individual was 
monitored. 


Statistical Analyses 

We used basic descriptive and multivariate statistics to examine 
feeder visitation parameters. These visitation parameters may be 
considered a proxy for food resource consumption, though our 
data do not allow for the exact relationship between consumption 
and visitation to be explored. Note for logistical reasons data 
collection periods for camera trap and PIT-tag data differed, 
hence in the analyses treatment of season differed for each. 
Specifically, because collection of camera trap data encompassed 
two winters, seasonal period (season ID) was used for analysis of 
this dataset to enable variation between winters to be examined 
(in the absence of sufficient data to include year as a fixed effect). 
Means are shown with their standard errors (SEs) (x + SEM), and 
the critical a level was 0.05 for all tests. 


Species Visitation and Co-occurrence (Camera Trap 
Data) 

To explore species associations at feeders (camera trap data) 
we initially used a probabilistic model approach (Veech, 
2013) to test for overall patterns of co-occurrence between 
species pairs (across all sites and seasons). This approach, 
implemented in the cooccur package (Griffith et al., 2015) 
in R 3.4.0 (R Core Team, 2017), uses presence/absence data 
to calculate an expected frequency of co-occurrence between 
species pairs if they were distributed independently of one 
another across sites, or observations in this case. The model 
then calculates the probability that the observed co-occurrence 
frequency is greater than the expected frequency (a positive 
co-occurrence association), less than the expected frequency 
(negative association), or random. Here the model tests the 
probability of co-occurrence at the level of the visitation event. 
For this analysis we used feeder-visiting species whose occupancy 
across study properties was >50% (Supplementary Table 2; 
determined via bird surveys conducted over the study duration; 
Galbraith et al., 2015). 

We then examined the composition of and variation in bird 
assemblages visiting feeding stations simultaneously using the 
PERMANOVA-4 add-on to PRIMER (Anderson et al., 2008). 
We calculated Bray-Curtis distances on fourth-root-transformed 
data as there were large differences in the baseline abundance 
of species. We then performed a non-metric multidimensional 
scaling (NMDS; Kruskal, 1964) on the distance-matrix which is 
displayed as a reduced space plot where each point represents 
the bird species composition at each feeding site and time point. 
To visualise species associations we superimposed vector lines, 
the length of which represents how much weight a species 
carried in determining the position of the points in the plot. 
The degree to which two lines are aligned with each other 
show the extent to which two species are correlated, with lines 


pointing in opposite directions meaning negative correlation. 
To explore the impact of seasonal period and food availability 
on species composition at feeding stations, a permutational 
analysis of variance (PERMANOVA) was used. Initially we fitted 
a PERMANOVA model with four factors: season ID (winter 2012, 
spring 2012, summer 2013, and winter 2013), food remaining 
(bread + seed remaining, bread remaining, seed remaining, no 
food remaining), property identity (ID), and observation day 
(1, 2, or 3). Observation day was found to be non-significant 
with minimal variance and was removed from the final model. 
Property ID was treated as a random factor while season ID and 
food availability were treated as fixed factors, and all factors were 
crossed. 

To investigate predictors of feeder use at the species 
level we used a mixed model approach, focusing on the six 
most frequent feeder-visiting species in the camera trap data 
(Supplementary Table 3). Generalized Linear Mixed Models 
(GLMMs) implemented in R were used, accounting for the 
repeated measures structure of the data. Two feeder-use 
parameters, abundance at feeders and number of daily visitation 
events (daily visitation rate), were modelled as the response 
variables for each species. Abundance models tested whether the 
number of conspecifics (individuals of the same species) during 
a visitation event was influenced by the abundance of other 
visiting species and/or varied among seasonal periods. Visitation 
rate models tested whether the number of daily visitation events 
for a focal species varied among seasonal periods. In abundance 
models, we included season ID, food remaining, the abundances 
of each of the other five focal species in the visitation event, 
minimum daily temperature (°C), and rainfall (mm) as fixed 
effects. In visitation rate models, we included season ID, duration 
food available (min day~!), minimum daily temperature, and 
rainfall as fixed effects. Temperature and rainfall data were 
obtained from the NIWA National Climate Database (for the 
Albany, North Shore, Auckland weather station; https://cliflo. 
niwa.co.nz/, accessed 15 November 2016). While avian foraging 
patterns typically vary with time of day, food was not consistently 
available throughout the day in our study due to depletion 
(reflecting the typical feeding practices of the public). As such, 
feeder visitation across the day was likely to reflect food levels 
remaining rather than daily patterns of foraging activity. We 
accounted for this in abundance models by including the “food 
remaining” parameter, with minutes after sunrise also included as 
a fixed effect to account for expected variation in foraging activity 
over the course of the day. Similarly, for visitation rate models 
we calculated the duration that food was available each day, and 
included it as a fixed effect to account for its likely influence on 
daily visitation rates. We included property ID in all models as a 
random effect, to account for variation between sites. In initial 
models we also included day length (this was correlated with 
season ID and so was removed), as well as observation day in the 
random effects structure (the effect was negligible so data were 
pooled across observation days for final models). 

Prior to model fitting we assessed the distribution of each 
response variable using the fitdistrplus package (Delignette- 
Muller and Dutang, 2015) in R; the best-fitting distributions were 
used in the corresponding models. We estimated the parameters 
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using negative binomial GLMMs, fitted using a Laplace 
approximation of maximum likelihood in the glmmADMB 
package (Skaug et al., 2015), in all cases [with the exception 
of the Eurasian blackbird (Turdus merula) abundance model 
where a Poisson error structure was the better fit]. For initial 
analyses, both ordinary and zero-inflated models were fitted. 
We used Pearson residual plots and Akaike’s Information 
Criterion (AIC) to compare model fits, retaining the best fitting 
model for interpretation of effects. Additionally, we checked 
for overdispersion by dividing the sum of squared Pearson 
residuals by the residual degrees of freedom and comparing 
this to a x? distribution (Venables and Ripley, 2002). Zero- 
inflation improved model fit for the house sparrow abundance 
model only. In initial model fitting, the blackbird, common 
starling (Sturnus vulgaris), and silvereye (Zosterops lateralis) 
abundance models all had large SEs for the “seed” level of the 
“food remaining” parameter, as did the silvereye abundance and 
visitation rate models for the “summer 2013” level of the “season 
ID” parameter. This was due to separation in the data due to cells 
with zero frequencies within the contingency table of response 
variable x food available (and x season for silvereye). To yield 
sensible parameter SEs, we added a dummy row to the dataset 
that added a small non-negative constant (a count of 1) to the 
cells with zero counts (Agresti, 2002; Jones et al., 2012) with 
other cells containing mean values. We then refitted these models 
and checked parameter SEs. For the final models we assessed the 
significance of whole model terms using likelihood ratio tests 
(LRTs) implemented with the “drop1” function in R. Post hoc 
pairwise comparisons among levels of season ID were conducted 
using the multcomp package (Hothorn et al., 2008). 


Individual Visitation (PIT-Tag Data) 
We used GLMMs implemented in R to explore feeder 
visitation patterns of individuals (PIT-tag data). Two feeder- 
use parameters, daily presence at feeders (binomial; present 
= 1, absent = 0) and daily visitation rate (count data; total 
reads per day active at feeder), were modelled as the response 
variables. Daily presence data were modelled with a binomial 
error structure (logit link) using the Ime4 package (Bates 
et al., 2014). Daily visitation rate data were modelled using 
a negative binomial error structure using the glmmADMB 
package, after assessing the distribution using the fitdistrplus 
package. For both response variables we fitted initial models 
that included species as a fixed effect and individual ID as a 
random effect, to confirm differences among species in feeder 
use. This also allowed us to estimate the relative contribution 
of among-individual variability (i.e., individual heterogeneity) 
and among-species variability to the overall variation in feeder 
visitation data. We calculated the marginal R? (proportion of 
variance explained by the fixed effects) and the conditional 
R? (proportion of variance explained by fixed and random 
effects combined) using the “r.squaredGLMM” function in 
the R package MuMIn (Barton, 2015), which implements the 
methods of Nakagawa and Schielzeth (2013). R? 


(CONDITIONAL) 
-Rm ARGINAL) gives the random effect component of the 
variance, RRANDOM) RRANDOM) is analogous to repeatability 


(Nakagawa and Schielzeth, 2010; Crates et al., 2016). Note, 
estimating the variance components of fixed vs. random effects 
for negative binomial models required refitting models using 
the “glmmPQL” function of the MASS package (which uses 
penalized quasilikelihood for parameter estimation; Venables 
and Ripley, 2002), as glmmADMB objects cannot be passed 
to “r.squaredGLMM.” Parameter estimates were comparable 
between “glmmPQL’ and “glmmadmb” fitted models (see 
Supplementary Tables 4, 5); coefficients of the latter are presented 
in the results for consistency. 

We then fitted separate GLMMs for each species for each 
response variable to examine the effects of season on feeder- 
use parameters. For both daily presence and daily visitation 
rate models, we included season (winter, spring, summer, and 
autumn), minimum daily temperature, rainfall, and day (number 
of days since feeding started) as fixed effects. Individual ID was 
included in the models as a random effect to account for repeated 
measures from the same individuals and to determine the 
contribution of individual heterogeneity to model variance. For 
all species except house sparrow, there were too few individuals at 
each property to adequately estimate among-property variation, 
thus data were pooled across properties. For house sparrow data, 
larger sample sizes at each property enabled initial models to be 
fitted with property ID (Feeder 1, 2, and 3) included as a fixed 
effect. As property ID did not contribute significantly to model 
fit for either response variable (see Supplementary Table 6), we 
pooled data across properties for final analyses as per the other 
species. Significance of whole model terms was assessed using 
LRTs. 


RESULTS 


Camera Trap Visitation Events 

A total of 3066 visitation events were captured by the camera 
traps over 132 trap days. We recorded 725 visitation events in 
winter 2012, 1006 in spring 2012, 581 in summer 2013, and 
754 in winter 2013. There were 723 visitation events made 
when bread and seed were both available, representing the time 
period immediately following food provision. There were fewest 
observations when only seed remained (n = 31); typically this 
food type was depleted fastest, leaving only bread remaining 
(n = 839 observations). The other visitation events occurred 
when there was no food remaining (n = 1473), although a small 
amount may have still been available on the table. 


Structure of Bird Assemblage at Feeding 


Stations 

Eleven bird species were recorded visiting the feeding stations 
(Supplementary Table 3). House sparrows and spotted doves 
were the most frequently observed species over all visitation 
events, present in 64.9 and 58.0% of recordings, respectively. 
Silvereyes were the only native species, present in only 4.6% of 
visitation events. Most visitation events had only a single species 
present (57.3%), with two species present in 33.5% of visitation 
events. Only 9.2% of visitation events had three or more species 
present. Mean total abundance per visitation event was 4.9 + 
1.53 birds, with a mean species richness of 1.5 + 0.01 species. 
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A maximum of five species were observed together on the feeding 
station in one visitation event. Mean conspecific group size 
(number of same-species individuals feeding simultaneously) 
was highest in house sparrow with 5.2 + 0.12 individuals per 
visitation event, and a maximum of 45 sparrows recorded in one 
visitation event. In contrast, Eurasian blackbirds, song thrushes 
(Turdus philomelos), and chaffinches (Fringilla coelebs) were 
typically observed without conspecifics present (Supplementary 
Table 3). 





Species Visitation Patterns 

House sparrows and spotted doves had the highest daily 
visitation rates (number of visitation events recorded by camera 
traps) to feeding stations (mean no. visitation events day! 
= 15.08 + 0.71 and 13.46 + 0.63, respectively), well above 
that for any other species (Supplementary Table 3, Figure 1A). 
Abundances of all six focal species (the most frequent feeder- 
visitors) at feeding stations were significantly influenced by 
food levels remaining (GLMM analyses; Table 1). In particular, 
abundances tended to be lowest when no food remained in 
comparison to when bread and seed remained (Wald-Z < 
—3.01, p < 0.003 in all cases). Furthermore, the speed at 
which food was depleted varied significantly among seasons 
(GLMM LRT: x? = 28.2, d.f. = 3, p < 0.001); food was 
available for longer in winter 2012 (mean min available day—! 
= 247 min + 27) and spring 2012 (290 min + 30) compared 
to summer 2013 (134 min + 17), when food was depleted 
fastest, and winter 2013 (170 min + 18). This duration of 
food availability significantly affected the number of daily 
visitation events by Eurasian blackbird and common myna, 
and marginally improved model fit for spotted dove and 
silvereye daily visitation rate models (Table 2). Nevertheless, after 
accounting for this variation in food availability, abundances 
at feeding stations and daily visitation rates for all species 
varied significantly among seasonal periods (except the spotted 
dove abundance model; Tables 1, 2 and Figures 1A,B). This 
seasonality was most striking for silvereyes, with no feeder 
visitations recorded for the summer period at all (Figure 1). 
House sparrow abundance at feeders was higher in winter 
2013 compared to winter 2012, whereas both daily visitation 
rate and abundance were lower in the second winter for 
Eurasian blackbirds (Figure 1, Supplementary Table 7). The 
remaining species showed no significant difference between 
winters for either visitation parameter (Figure 1, Supplementary 
Table 7). 








Species Associations 

A number of significant associations in overall species co- 
occurrence patterns at feeding stations were detected with the 
probabilistic modelling (Figure 2). Of the 36 species pairs, three 
were positive (8.3%) with the two species co-occurring at feeders 
significantly more frequently than expected, 10 were negative 
(27.8%) with species co-occurring significantly less frequently 
than expected, and 23 were random (63.9%). House sparrow had 
the highest number of significant associations with other species 
(n = 6), and the highest proportion of positive associations 
(37.5% of pairings). Eurasian blackbirds had five significant 


associations with other species, while spotted doves had four, 
with all of these being negative. 

Analysis at the community level also indicated there were 
patterns in the assemblages of birds feeding concurrently at the 
feeding stations. Species composition was dominated by three 
species, house sparrows, spotted doves, and Eurasian blackbirds, 
with strong negative correlations between the abundances of 
these species (Figure 3); blackbirds were present at feeding 
stations when there were less spotted doves and house sparrows. 
Species composition at feeding stations varied significantly with 
seasonal period and food availability, and was influenced by 
property ID (PERMANOVA: F = 2.38, d.f. = 54, P = 0.001). 

GLMM analyses provided further evidence of species 
associations, with the abundance of all focal species during 
visitation events influenced by the abundance of at least one other 
top feeder-visitor (Table 1). Significant effects were generally 
negative, with the exception of the house sparrow and silvereye 
models. With all other factors held constant, house sparrow 
abundance increased with common starling, silvereye, and 
common myna abundances, and silvereye abundance increased 
with house sparrow abundance. Effect sizes for season ID and 
food remaining were typically larger than for co-occurring 
species predictors though. 


Individual Visitation Patterns 

We PIT-tagged a total of 110 individuals from five feeder- 
visiting species (Supplementary Table 1). Numbers tagged were 
equivalent at each feeding station (Feeder 1: n = 37; Feeder 2: n = 
37; Feeder 3: n = 36). The redetection rates varied among species 
(Supplementary Table 1), with 70 individuals (63.6%) overall 
redetected by RFID readers at feeding stations on at least one 
occasion after initial capture. In total 83265 reads were recorded 
(with duplicate reads removed) for these 70 individuals. 

There were obvious differences among individuals in the 
consistency of feeder visitation over time, with some individuals 
returning almost daily while others only visited sporadically 
(Supplementary Figure 3). Daily presence at feeding stations 
differed significantly among species (i.e., the proportion of days 
birds returned to the feeder out of the total days monitored; 
GLMM LRT: x? = 18.9, d.f. = 4, P < 0.001; Figure 4; 
Supplementary Table 4). The variation among individuals 
explained a greater proportion of the variation in daily 
presence data than among-species variation (48.4 vs. 14.8%, 
respectively). Considering each species separately (except for 
common starling for which there were too few individuals for 
valid comparisons), GLMM analyses indicated that the daily 
likelihood of individuals visiting feeders varied seasonally for all 
species except house sparrow (Table 3). House sparrows showed 
the greatest variability among individuals, with individual 
heterogeneity accounting for 63.2% of the variance in the model 
[Eurasian blackbirds 52.1%, spotted doves 38.4%, Barbary doves 
(Streptopelia roseogrisea) 14.3%; Table 3]. 

Mean daily visitation rates varied significantly between species 
(GLMM LRT: x? = 17.6, d.f. = 4, P < 0.01; Supplementary 
Table 4). However, heterogeneity among individuals explained 
a greater proportion of the variance in daily visitation rate 
than among-species variation (19.6 vs. 4.9%, respectively). Three 
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FIGURE 1 | Patterns of feeder visitation of the top six avian visitors captured in camera traps at experimental feeding stations at urban study properties in northern 
Auckland, New Zealand. Mean count data (+ SEM) are presented for two visitation parameters scored from camera trap data: (A) daily visitation rate (number of daily 
visitation events recorded); and (B) abundance (number of individuals in a visitation event). In each panel, means not sharing the same letter are significantly different in 
pairwise comparisons (Tukey’s HSD; p < 0.05; Supplementary Table 7). Sampling season ID: W1, winter 2012; SP, spring 2012; SU, summer 2013; W2, winter 2013. 
Note y-axis scale varies with species. 











house sparrow individuals had the highest visitation rates per season for all species (Table 4). Individual heterogeneity in daily 
active day, with means of 32.7 (+ 14.5) to 42.8 (4 8.2) visitation rate was highest for house sparrows and blackbirds, 
reads per day (Figure 5). Analysed separately (again except for accounting for 31.1 and 28.4% of explained variance in the 
common starling), daily visitation rates varied significantly with models, respectively (Table 4). 
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TABLE 1 | Results of Generalized Linear Mixed Models (GLMMs) examining factors affecting the abundance of garden birds captured in camera traps at experimental 
feeding stations in northern Auckland, New Zealand. 























Explanatory variables House sparrow Spotted dove Eurasian blackbird Common starling Silvereye Common myna 
Intercept 1.806 (0.143) 0.957 (0.188) —1.179 (0.275 —1.713 (0.433) —1.439 (0.644) —1.811 (0.449) 
Season ID (reference: Winter 2012) a ° S ai ue = 

Spring 2012 —0.255 (0.073) 0.064 (0.063) 0.001 (0.143) —0.116 (0.257) —2.014 (0.300) —0.278 (0.330) 

Summer 2013 0.726 (0.103) —0.062 (0.099) —0.579 (0.258 —1.829 (0.430) —3.975 (1.084) —1.993 (0.588) 

Winter 2013 0.262 (0.074) 0.103 (0.067) —0.805 (0.216 —0.415 (0.263) 0.375 (0.277) —0.383 (0.346) 
Food remaining (reference: Bread and Seed) Ta — ii 7 ue ii 

Seed only 0.668 (0.177) —0.712 (0.262) —0.198 (1.034 —1.402 (1.094) —1.768 (1.158) —1.081 (0.909) 

Bread only —0.340 (0.064) —0.909 (0.053) 0.038 (0.162) 0.216 (0.207) —0.489 (0.211) —0.604 (0.306) 

None —1.094 (0.067) —1.223 (0.059) —0.721 (0.175 —0.807 (0.268) —1.776 (0.289) —2.397 (0.450) 
House sparrow abundance - —0.008 (0.005) * = —0.064 (0.018)*** 0.073 (0.018) 0.040 (0.018)* 0.038 (0.025) 
Spotted dove abundance —0.018 (0.015) - —0.590 (0.076)*** —0.060 (0.049) —0.292 (0.071)*** — —0.088 (0.063) 
Eurasian blackbird abundance —0.353 (0.072)*** —0.835 (0.081)*** - —0.458 (0.295) —0.821 (0.331)** —2.135 (0.749)*** 
Common starling abundance 0.101 (0.050)* 0.020 (0.048) —0.373 (0.206)* - —0.740 (0.410)* —0.366 (0.276) 
Silvereye abundance 0.117 (0.048)* —0.420 (0.063)""* —0.615 (0.228)*** —0.677 (0.402)* - —0.324 (0.253) 
Common myna abundance 0.179 (0.078)* —0.074 (0.056) —1.345 (0.531)*** —0.488 (0.304) ° —0.055 (0.287) - 
Time of day (min after sunrise) —0.001 (0.000)*** —0.001 (0.000)*** 0.0003 (0.000) —0,.004 (0.001)*** = —0.001 (0.001) —0.003 (0.001)** 
Minimum daily temperature (°C)* 0.001 (0.009) 0.007 (0.008) —0.009 (0.019) 0.0004 (0.033) —0.082 (0.034)* 0.059 (0.046) 
Rainfall (mm)* 0.009 (0.006) 0.012 (0.004)** 0.034 (0.017) ° —0.060 (0.035) ° —0.003 (0.015) —0.004 (0.030) 
Property ID (random effect)t 0.132 (0.364) 0.341 (0.583) 0.410 (0.641) 1.119 (1.058) 3.136 (1.771) 0.297 (0.545) 








Parameter estimates and their standard error (+ SE) are presented for each model term at the reference levels stated. Whole effects were tested with likelihood ratio tests (LRTs); 
significant chi-square test statistics from LRTs are indicated with: ep < 0.10, *p < 0.05, “p < 0.01, and ***, p < 0.001. 

Models were fitted with negative binomial error structures, except for house sparrow (zero-inflated negative binomial) and Eurasian blackbird (Poisson). 

* Data from NIWA National Climate Database. 

t Variance (standard deviation) presented. 

Models were fitted for the six most frequent feeder-visiting species (presented here in descending order of frequency). n observations = 3066. 


TABLE 2 | Results of Generalized Linear Mixed Models (GLMMs) testing whether the number of visitation events per day captured by camera traps varied among seasons. 











Explanatory variables House sparrow Spotted dove Eurasian blackbird Common starling Silvereye Common myna 
Intercept 2.512 (0.161) 2.204 (0.146 0.595 (0.405) 0.387 (0.481) —0.065 (0.666) 0.239 (0.393) 
Season ID (reference: Winter 2012) “= = j ii "o i 

Spring 2012 0.338 (0.088) 0.425 (0.108) 0.390 (0.278) 0.206 (0.335) —1.748 (0.413) 0.115 (0.319) 

Summer 2013 0.290 (0.130) 0.043 (0.161 —0.526 (0.457) —1.349 (0.528) —3.585 (1.167) —1.426 (0.564) 

Winter 2013 —0.130 (0.092) 0.023 (0.114) —1.119 (0.331) —0.520 (0.336) 0.058 (0.421) —0.734 (0.356) 
Duration food available (min day~') —0.0002 (0.0002) 0.0005 (0.0003) ° 0.002 (0.0008)** —0.0008 (0.0009) 0.002 (0.001) ° —0.002 (0.001)* 
Minimum daily temperature (°C)* 0.0003 (0.011) 0.012 (0.014) —0.028 (0.037) —0.002 (0.044) —0.033 (0.050) 0.014 (0.045) 
Rainfall (mm)* 0.005 (0.007) —0.001 (0.008) 0.018 (0.028) —0.057 (0.041) —0.008 (0.022) 0.010 (0.027) 
Property ID (random effect)t 0.189 (0.434) 0.077 (0.279 0.640 (0.800) 1.275 (1.129) 2.341 (1.530) 0.098 (0.027) 











Parameter estimates and their standard error (+ SE) are presented for each model term at the reference levels stated. Whole effects were tested with likelihood ratio tests (LRTs); 
significant chi-square test statistics from LRTs are indicated with: ep < 0.10, *p < 0.05, **p < 0.01, and ***p < 0.001. Models were fitted with negative binomial error structures. 

# Data from NIWA National Climate Database. 

t Variance (standard deviation) presented. 

Models were fitted for the six most frequent feeder-visiting species (presented here in descending order of frequency). n observations = 132. 


DISCUSSION 


Summary of Findings 

Our results demonstrate considerable among- and within-species 
variation in the use of supplementary food resources by birds 
in situ in urban habitats, confirming that feeder-visiting birds 
do not exploit supplementary food resources equally. The 
species using feeding stations in our study were predominantly 


introduced; this was predicted because of the grain-based feeding 
regime used by householders (in our study and more generally; 
Galbraith et al., 2014), and the division in dietary preferences 
of native vs. introduced birds in our study area (Galbraith 
et al., 2015). Two species in particular dominated resource use: 
house sparrow and spotted doves. Our study corroborates the 
findings of two other recent experimental studies, one focusing 
on great tits (Parus major) and blue tits (Cyanistes caeruleus) 





Frontiers in Ecology and Evolution | www.frontiersin.org 8 August 2017 | Volume 5 | Article 81 


Galbraith et al. 


Species-Specific Dominance at Bird Feeders 











FIGURE 2 | Pairwise associations between species at experimental feeding 
stations at urban study properties in northern Auckland, New Zealand, 
determined through probabilistic modelling of species co-occurrences (Veech, 
2013). Species that co-occur more frequently than by chance (a threshold = 
0.05) are considered to have a positive association (in blue), those that 
co-occur less frequently than by chance are considered to have a negative 
association (in red), with random associations shown in grey. Note chaffinch 
and greenfinch pairings are not shown; no significant associations were found 
for these species. Species masses are means from birds captured in this study. 
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FIGURE 3 | NMDS ordination of avian community composition at 
experimental feeding stations at urban study properties in northern Auckland, 
New Zealand. Each point (blue circles) represents the avian assemblage 
recorded in a single visitation event. The overlaid species vectors (red lines) 
illustrate that three bird species, Eurasian blackbirds (BLKB), house sparrows 
(SPRW), and spotted doves (SPDV) dominate community composition (as they 
have the longest vectors), and that there are strong negative correlations 
between the abundances of these species. BADV, Barbary dove; CHFN, 
chaffinch; GRFN, greenfinch; MYNA, common myna; RPIG, rock pigeon; 
SEYE, silvereye; STRL, common starling; THSH, song thrush. 











in woodland habitat in Oxford, UK, during winter (Crates 
et al., 2016) and the other examining visitation patterns of six 
species at a university campus in Cornwall, UK (Jack, 2016). 


100 


N 
a 


Daily presence at feeding station 
(% of days monitored) 
N a 
a oO 


Barbary dove ee eb o © 


House sparrow 
Spotted dove 


Common starling 
Eurasian blackbird 


Species 


FIGURE 4 | Consistency of feeder visitation of 70 individually PIT-tagged birds 
at experimental feeding stations in northern Auckland, New Zealand. Points 
represent the percentage of days an individual was active at a feeding station 
out of all days monitored for that individual. Horizontal lines represent the 
group means for each species. 











These studies both reported individual and interspecific variation 
in use of supplementary food at experimental feeders over a 
comparatively shorter timeframe (c. 3 months). The findings 
from our longer-term study also demonstrate that season has an 
important influence on a number of visitation parameters both at 
the individual and species levels. Relatedly, seasonal differences 
in movement of birds among networks of urban feeders have 
been found in great tits and blue tits (Cox et al., 2016). We 
acknowledge, though, that in our study there are limitations to 
interpreting the results from those PIT-tagged species with small 
sample sizes. Nevertheless, significant seasonal variation in feeder 
use was evident across most species and visitation parameters we 
assessed. 


Visitation Patterns 

The interspecific differences in feeder use observed in our 
study reinforce the concern that typical feeding practices in 
New Zealand do not benefit native species, but instead support 
introduced species (Galbraith et al., 2015). Feeding grain-based 
food types (e.g., bread, seed) is the prevalent practice in New 
Zealand (Galbraith et al., 2014), yet most native garden birds 
here are frugivores, nectarivores, or insectivores (Heather and 
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TABLE 3 | Results of Generalized Linear Mixed Models (GLMMs) testing the effect of season on daily visitation (presence/absence) of individually PIT-tagged birds to 
experimental feeding stations. 














House sparrow Spotted dove Eurasian blackbird Barbary dove 

n observations 13367 1453 2088 2479 
n individuals 45 7 7 7 
Intercept 0.435 (0.423) 4.969 (0.834) 0.168 (0.999) —0.836 (0.841) 
Season (reference: Winter) = e d 

Spring —0.058 (0.109) —0.156 (0.259) 0.675 (0.280) 1.599 (0.488) 

Summer 0.143 (0.112) —1.050 (0.337) 1.043 (0.280) 1.403 (0.451) 

Autumn 0.090 (0.082) —0.994 (0.218) 0.345 (0.213) 0.165 (0.424) 
Minimum daily temperature (°C)* 0.003 (0.009) —0.031 (0.026) 0.028 (0.022) —0.004 (0.026) 
Rainfall (mm)* —0.009 (0.004)* 0.001 (0.009) —0.006 (0.009) —0.026 (0.018) 
Day (numerical from start of feeding) —0.004 (0.000)*** —0.007 (0.001)*** —0.007 (0.001)*** —0,.008 (0.001)*** 
Individual ID (random effect)t 5.992 (2.448) 2.388 (1.545) 4.767 (2.183) 0.906 (0.952) 
Marginal R2 0.021 0.088 0.119 0.337 
Conditional R2 0.653 0.471 0.640 0.480 








Parameter estimates and their standard error (+ SE) are presented for each model term at the reference levels stated. Whole effects were tested with likelihood ratio tests (LRTs); 
significant chi-square test statistics from LRTs are indicated with: ep < 0.10, *p < 0.05, **p < 0.01, and ***p < 0.001. Models were fitted with binomial error structures. 

# Data from NIWA National Climate Database. 

t Variance (standard deviation) presented. 

Species ordered by frequency of occurrence from Supplementary Table 3. 
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FIGURE 5 | Daily rates of feeder visitation (number of visitation events per day actively using the feeder) of 70 individually PIT-tagged birds at experimental feeding 
stations in northern Auckland, New Zealand. Data are raw means + SEM. 











Robertson, 1996). Thus, unsurprisingly, native species were experiment, in which house sparrow and spotted dove showed 
almost entirely absent from our experimental feeding stations, dramatic increases in abundance at feeding locales in response 
with the exception of the silvereye which primarily visited in to the feeding treatment, contributing to a shift in community 
winter. Of the 11 recorded feeder-visiting species, two largely | composition (Galbraith et al., 2015). This is concerning for 
granivorous introduced species—house sparrow and spotted urban ecosystems in New Zealand, where introduced birds 
dove—dominated the majority of visitation events. They were are already the dominant component of avifaunal assemblages, 
also the two most abundant species at feeders, both in terms and increasingly relevant globally, as supplementary feeding is 
of mean and maximum conspecific group size using the implicated as a contributing factor in the spread of other invasive 
feeder simultaneously. These results reflect the findings of bird bird species, such as monk parakeets (Myiopsitta monachus) in 
counts conducted at study properties over the duration of the the northern USA (Davis et al., 2014) and ring-necked parakeets 
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TABLE 4 | Results of Generalized Linear Mixed Models (GLMMs) testing the effect of season on daily visitation rate (no. visits per active day) of individually PIT-tagged 


birds at experimental feeding stations. 








House sparrow Spotted dove Eurasian blackbird Barbary dove 

n observations 4695 990 787 316 
n individuals 45 IG 7 7 
Intercept 3.423 (0.144) 4.831 (0.205) 4.145 (0.395) 5.735 (0.812) 
Season (reference: Winter) — =~ m iii 

Spring —0.052 (0.054) 0.242 (0.086) —0.221 (0.119) —0.740 (0.558) 

Summer —0.469 (0.056) —0.403 (0.117) —1.052 (0.124) —1.746 (0.508) 

Autumn —0.054 (0.039) —0.055 (0.072) —0.913 (0.105) —1.049 (0.511) 
Minimum daily temperature (°C)# —0.012 (0.004)** —0.036 (0.008)*** —0.042 (0.010) —0.000 (0.025) 
Rainfall (mm)* 0.001 (0.002) 0.001 (0.003) 0.008 (0.005) —0.022 (0.018) 
Day (numerical from start of feeding) —0.003 (0.000)*** —0.005 (0.000)*** —0.005 (0.000)*** —0.011 (0.001)*** 
Individual ID (random effect)t 0.428 (0.654) 0.067 (0.259) 0.672 (0.820) 0.037 (0.193) 
Marginal R2 0.065 0.197 0.226 0.505 
Conditional R2 0.376 0.247 0.510 0.522 








Parameter estimates and their standard error (+ SE) are presented for each model term at the reference levels stated. Whole effects were tested with likelihood ratio tests (LRTs); 
significant chi-square test statistics from LRTs are indicated with: ep < 0.10, *p < 0.05, “p < 0.01, and **p < 0.007. 


Models were fitted with negative binomial error structures. 

# Data from NIWA National Climate Database. 

t Variance (standard deviation) presented. 

Species ordered by frequency of occurrence from Supplementary Table 3. 


(Psittacula krameri) in the UK and mainland Europe (Strubbe 
and Matthysen, 2007; Peck et al., 2014; Orros and Fellowes, 
2015b). This concern is not limited to introduced birds, however, 
as native species may also become numerically dominant as a 
result of food subsidies and displace other birds. 

The prevalence of house sparrows and spotted doves in 
visitation events and their distinct positive numerical responses 
to the feeding treatment (Galbraith et al., 2015) suggests 
both species were successfully competing for access to the 
food resource. However, there were differences between house 
sparrows and spotted doves in visitation patterns which may 
be indicative of distinct, but equally successful, competitive 
behaviours. Spotted dove feeder visitors appear to consist of a 
core-group of individuals that are highly consistent in their feeder 
use (i.e., show high fidelity to the resource once they discover it), 
whereas house sparrow feeder visitors are highly numerous but 
also highly variable in their feeder use. 

Our results provide multiple lines of evidence (both camera 
trap and PIT-tag data), that seasonality has an important 
effect on feeder visitation patterns. Unsurprisingly, the level of 
food remaining had a significant impact on feeder visitation 
parameters—when there was more food left, bird abundance at 
feeders was higher. Nevertheless, after accounting for the effect 
of food availability, species abundance and species daily visitation 
rate varied significantly with seasonal period for most focal 
species. We also found that the daily presence and visitation rates 
of individuals varied significantly among seasons for all focal 
species (though we could not account for food availability here). 
There were, however, two exceptions: first, the number of spotted 
doves feeding at the same time was consistent among seasonal 
periods. This may well reflect the maximum capacity of the 
feeding station structure itself—ultimately the number of birds, 
particularly large-bodied species, such as spotted doves, feeding 


simultaneously will be restricted by physical space at feeders. 
Second, no seasonal effect on the likelihood of individuals being 
present each day was evident for house sparrow—individual 
heterogeneity was more influential compared to season for this 
species. 

The physiological demands on birds vary seasonally with 
changes in environmental conditions, natural food availability, 
and behavioural activity, thereby influencing the demand for 
supplementary food. Here our proxy for demand was the speed of 
food depletion, which indicated that demand for supplementary 
food resources did fluctuate seasonally in our study. Demand 
was highest in summer rather than winter, which is contrary 
to the common public perception that birds should be fed in 
the winter when their struggle for survival is greatest (Jones, 
2017). However, in our study region winter conditions are fairly 
mild, so the physiological demand is likely less compared to 
many regions. Instead, the energetic costs of breeding activities 
(e.g., increased energy required for reproduction, and/or reduced 
foraging time due to incubation and chick provisioning activities; 
Martin, 1987) may drive the increased demand for readily 
available supplementary food in summer. Additionally, it is 
crucial to note that our study period encompassed a pronounced, 
New-Zealand-wide drought event over the summer of 2012- 
2013 (Porteous and Mullan, 2013). This event would probably 
have had significant impacts on natural food availability, and 
may well have contributed to food demand. The differences 
observed between winters may also reflect follow-on impacts 
of the drought event, if, for instance, natural food resource 
availability had not recovered. However, visitation patterns 
of the only native feeder visitor provide a notable contrast 
to the overall seasonal pattern of demand. Silvereyes were 
scarce at feeders in spring and entirely absent in summer, 
suggesting that natural food resources exploited by this species 
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were more abundant (or more profitable) during these seasons 
and that supplementary food was unnecessary to meet energy 
requirements. 


Species Associations 

We found evidence that access to supplementary food was 
influenced by species associations at feeding stations. Feeder 
activity typically peaked in the morning following food 
provisioning. Although multiple species were typically present 
in the vicinity of feeders, single-species visitation events were 
more common than those with multiple species at the feeder 
simultaneously. It is likely, then, that some species excluded 
others from access to feeders and some avoided using feeders 
when heterospecifics were present. Multi-species visitation 
events, however, still comprised a substantial proportion of our 
observations (42.7% of visitation events). 

A number of significant interspecific associations were 
identified from the probabilistic modelling of the likelihood of 
co-occurrence between species pairs, which generally agreed with 
findings from GLMM analyses of predictors of species abundance 
at feeders. Body size frequently determines, or contributes to, the 
development of dominance hierarchies, outcomes of agonistic 
interactions, and subsequent resource access (Forrester, 1991; 
Wojczulanis-Jakubas et al., 2015; Miller et al., 2017). However, 
there were no obvious patterns in the species co-occurrence 
matrix to indicate that negative associations were more common 
between species closest in body size, or that positive associations 
were more common between species of greater size disparity. 
Thus, potentially, behavioural differences among species (i.e., 
in foraging strategies, territoriality, and sociality), rather than 
morphological differences alone, are driving co-occurrence 
patterns in this system. 

The known ecology of the focal species provides some insight 
here. For instance, house sparrow are a highly gregarious species, 
typically foraging, roosting, and nesting in conspecific flocks, 
and frequently foraging with other species (Higgins et al., 
2006a). Consistent with this apparent tolerance of foraging 
with heterospecifics, house sparrow had the highest number 
of positive associations with other feeder-visitors. In contrast, 
Eurasian blackbirds are highly territorial, typically forage singly 
or in pairs, with agonistic interspecific interactions occurring 
where they have been observed foraging with other species 
(Higgins et al., 2006a). This was reflected in our observed mean 
conspecific group size for blackbirds at feeders (1.04 birds) 
and the predominantly negative associations we found between 
blackbirds and other feeder-visiting species. Thus, blackbirds 
were apparently intolerant of using feeders in the presence of 
other species or could not gain access to them. They may 
have avoided using the feeder when other birds were present, 
actively excluded other species, been excluded by other species, 
or a combination of these. Detailed behavioural analysis of 
camera trap videos would help elucidate the direction of these 
interspecific interactions. 


Underlying Foraging Mechanisms 
Although we did not explicitly investigate which competitive 
mechanisms controlled feeder access, it seems probable that 


both interference and exploitative competition were at play. 
Agonistic interactions at feeding stations were regularly observed 
in this study (between species but more frequently within species; 
JAG pers. obs.) and are frequent among feeder-visiting birds in 
other systems (e.g., Wojczulanis-Jakubas et al., 2015), indicating 
that interference competition certainly influences feeder access 
to some degree. However, in our study system exploitative 
competition is also likely to be a critical mechanism. For 
example, house sparrows and spotted doves typically arrived at 
feeding stations in the greatest numbers and consumed food 
quickly, which are important indicators of exploitative ability 
(Holway, 1999; Bertelsmeier et al., 2015). Additionally, house 
sparrows forage in localised areas, typically moving through 
slowly (Higgins et al., 2006b); spending longer in a particular 
patch adds to efficiency in resource exploitation (Holway, 1999; 
Shochat et al., 2004). The exploitative ability of the spotted 
dove may also be enhanced by the reduced food handling time 
and food intake rates associated with consuming seeds whole 
(Shochat et al., 2004). 


Individual Variability 

Individual variation or “specialisation” is widespread in many 
taxa, yet conspecific individuals are treated as ecologically 
equivalent in the majority of studies concerning resource use 
(Bolnick et al., 2003; Araújo et al., 2011). Here we found 
convincing evidence that, within species, individuals varied 
considerably in their feeder use. Individual heterogeneity was 
highest in the house sparrow, in terms of both feeder visitation 
rates per day and likelihood of daily presence at feeders, with 
these individual differences highly repeatable. This high level of 
individual variation may result from a number of behavioural 
mechanisms. For instance, it may be indicative of an intraspecific 
dominance hierarchy dictated by the competitive ability of 
individuals (Richner, 1989; Rat et al., 2015). Alternatively, 
individual variation may reflect behavioural flexibility, a trait 
which has been linked to successful urban exploiter and 
invasive species (Sol et al., 2002; Peck et al., 2014), especially 
in regard to flexible foraging strategies. Distinctive individual 
foraging strategies have been found, for example, within white 
storks (Ciconia ciconia) in southwestern Spain which show 
both consistent (i.e., specialist) and flexible (i.e., generalist) use 
of landfills vs. rice fields (Sanz-Aguilar et al., 2015). These 
differences in foraging strategy within a species may arise 
from phenotypic, physiological, or competitive discrepancies 
among individuals (Araújo et al., 2011; Sanz-Aguilar et al., 
2015). Individual-level variation in feeder use may also be 
indicative of differences in personality (Aplin et al., 2014). 
Among Eurasian blackbird individuals the high variation in 
feeder use may reflect territorial behaviour. Expectedly, territory 
owners would use feeders more frequently and with greater 
consistency, whereas territory intruders would face a greater 
challenge accessing feeders due to territorial defence behaviour of 
the territory owner, and so be less frequent and/or less consistent 
feeder visitors. Studies in territorial calliope hummingbirds 
(Selasphorus calliope; Tamm, 1985) and great tits (Ydenberg, 
1984) have found that individual birds that have been provided 
with supplementary food show increased display and/or defence 
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behaviours, thus making it more difficult for others to access the 
resource. 

Individual heterogeneity in feeder use has important 
implications at both the individual-and population-level. Among 
other impacts, individuals that heavily rely on supplementary 
food may have altered body condition, reproductive success, 
home range size, and/or survival likelihood (Annett and 
Pierotti, 1999; Oro et al., 2013). The collective impacts of 
supplementary feeding on individuals determines the overall 
effect on population dynamics, including changes in population 
size, distribution, and migratory patterns (Robb et al., 2008; 
Oro et al., 2013; Amrhein, 2014), and will depend on both 
the proportion of the population using feeders and individual 
heterogeneity in feeder use, particularly where feeder users are 
a nonrandom subset of the population (Sanz-Aguilar et al., 
2015). Furthermore, individuals which follow a high-use, 
consistent foraging strategy at feeders are more likely to pose a 
disease risk than sporadic users, as exposure to pathogens and 
parasites increases as feeder visitation increases (Adelman et al., 
2015). 


CONCLUSIONS 


Our study highlights that individual and species-specific 
differences in feeder use are present within feeder-visiting bird 
communities, importantly demonstrating this across seasons 
within an urban system. These intraspecific and interspecific 
asymmetries support the likelihood of competitive interactions 
operating to regulate access to food, and suggest that the 
effects of supplementary feeding are unlikely to be equivalent 
across all birds within communities of feeder visitors. In 
New Zealand resource dominance by introduced species is 
particularly important, with negative outcomes for native species 
conservation in cities possible. Individual differences in feeder 
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